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Abstract. Crystallography has proven a rich source of ideas over several centuries. Among 
CO I the many ways of looking at space groups, N. David Mermin has pioneered the Fourier- 

■ space approach. Recently, we have supplemented this approach with methods borrowed 
I from algebraic topology. We now show what topology, which studies global properties of 

manifolds, has to do with crystallography. No mathematics is assumed beyond what the 
^ I typical physics or crystallography student will have seen of group theory; in particular, the 

■ reader need not have any prior exposure to topology or to cohomology of groups. 
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' 1. Introduction 

^ . Reviewing the Fourier-space formulation of crystallography, David Mermin wrote in 1992 

\^ [ that 

^ ! More than one person has told me that what I am calculating here are coho- 
rt i mology groups. I have found this information less valuable than M. Jourdain 
^ ■ found the news that he was speaking prose, but am too ignorant to state with 
^ • confidence that this is not a useful point of view. [22] 

^ ■ Applying cohomology, we have proven certain theorems 1^ 1^1 IE] in niore generality than 

. they were previously known, including to cases with 46-fold [21] and other exotic rotational 

Ch ! symmetries. Our purpose here, using minimal jargon and mathematical machinery, is to 

2 ' explain just what the homological point of view is. Since this point of view represents, in 

^ ■ some historical sense, a marriage between group theory and topology, and since most readers 

• 1— I . will already be familiar with the crystallographic implications of the former, we will display 

^ I some easily appreciated topological analogues and offer a complete example space-group 

d ■ calculation, without shortcuts, employing some of the new tools. 



Most concisely, topology is concerned with numbers that remain invariant under transfor- 
mations; for example, the number of holes in a doughnut does not change under rotations, 
continuous stretching, or twisting. In physics, electric charge does not change under a gauge 
transformation, field theories routinely refer to scalars as "charges," and the Kosterlitz- 
Thouless treatment of two-dimensional phase transitions has given rise to the description of 
vorticity as "topological charge" In the traditional, direct-space, formulation of crystal- 
lography, the space group of a structure does not depend on the choice of origin ("setting"), 
despite the fact that not all point-group elements need pass through the same point in a 
unit cell. There is a long history of the application of group cohomology to describe this 
invariance [21 13 ISH [131 dH • 

Similarly, Fourier-space crystallography, as first formulated by Bienenstock and Ewald 
[1] and developed by Mermin and collaborators jSHl I2H1 1231 1221 1211 IHl UH] , admits certain 
quantities invariant under a simpler analogue (to be defined below) of the electromagnetic 
gauge transformation. In the special case of a periodic crystal, this gauge transformation 



tWith apologies to N.D. Mermin, [HEOj- 
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can be described by a translation for each point-group symmetry. In the initial work on 
Fourier-space crystallography, tracking down these invariants was incidental to the main 
task of classifying space groups; most of the time, the invariant of a non-symmorphic space 
group corresponds to a necessary extinction in diffraction (something that clearly should not 
depend on the arbitrary choice of a gauge). Mermin first noticed that the invariants might 
not always be so simple in the very article alluding to the Bourgeois Gentleman.^ There, 
he noted that of the 157 non-symmorphic periodic space groups in three dimensions, and 
additionally of an important infinite class of quasiperiodic space groups that includes all 
the quasicrystals so far discovered, only two (both among the 157) had invariants that did 
not correspond to extinctions in diffraction.^ Later, Mermin and Konig |^ showed 
that the new kind of invariant had a different physical interpretation, that of electronic 
"band sticking." The Konig-Mermin conditions on the invariant for band sticking were 
too restrictive; we have generalized them |9j and have also found a third type of invariant 
corresponding to neither extinctions nor Konig-Mermin band sticking |H]. 

1.1. Fourier-Space Crystallography. In order to establish notation, we briefly review 
the Rokhsar-Wright-Mermin formulation of crystallography in Fourier space; see reference 
jn] or JH] for a more careful development. We begin with the (reciprocal) lattice L, which 
consists of all integral linear combinations of a finite number of generating vectors.^ The 
Fourier transform of a density (mass, electronic, etc.) has support on this lattice; we will 
refer to the Fourier transform as the density, since the direct-space density plays no role in 
the theory. Two densities, p(k) and p'(k), are indistinguishable if all their n-body correlation 
functions are the same; this is equivalent to the condition 

(1.1) p'(k) = e2-^>^(^)p(k) , 

where the gauge function x(k), defined modulo unity, is linear on the lattice. In case a point 
operation g takes p to p', we write the special gauge function associated with g as the phase 
function the point group G consists of all g such that 

(1.2) p(^k) = e2"'*«Wp(k) • 

The phase function is the central object of the Fourier formulation, and it is subject to a 
restriction and to an identification. As a consequence of the associativity of the group action 
on the lattice, the phase function must satisfy the group- compatibility condition for g,h & G, 

(1.3) <l>,^(k) = $g(/ik) + <l>;,(k) , 

^Two of us were among those who suggested to David Mermin that Fourier-space crystallography could 
be described in the language of group cohomology; Andre LeClair was another jl7| . As far as we can 
tell, the existence of this connection between Fourier-space crystallography and group cohomology was first 
mentioned in print in Mermin's Moliere reference; in the same year, Piunikhin alluded briefly to the same 
thing |25j . While it is perhaps interesting that results first formulated in one language can be recast in 
another, the correspondence becomes useful when it leads to new results. 

^The absence of systematic extinctions in these two space groups, /2i2i2i and /2i3, was well known. 
What was new was the description of an algebraic invariant that "detects" that these two groups are non- 
symmorphic. 

'^The lattice of a periodic structure is discrete and generated by three vectors in three dimensions, while 
that of a quasiperiodic structure is not discrete and has more generators than dimensions. In Fourier-space 
crystallography, this is the only theoretical distinction between these two cases. Since we are not interested 
in the direct-space lattice (if it even exists), we shall use the term "lattice" for the the reciprocal lattice. 
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while two phase functions, $ and describe indistinguishable densities if they are related 
by a gauge function, x, 

(1.4) $;(k) = $,(k) + x(^7k-k) , 

and so are identified. The lattice L, the action of the point group G, and the phase function <l> 
determine the space group type of the crystal. Homology theory provides a convenient way 
of calculating all such functions, subject to the restriction (jl.Hj) and the identification ()1.4|1 . 

2. Invariants 

We continue the review of Fourier-space crystallography, stressing the role played by gauge- 
invariant values of combinations of phase functions. We then switch context to elementary 
topology, where we define homology of loops drawn on two-dimensional surfaces in an in- 
tuitive way. Cohomology is introduced through a familiar example from vector calculus. 
Finally, we connect topology to crystallography by way of homology and cohomology of 
groups. 

2.1. Invariants in Fourier-space Crystallography. Evidently, gauge-equivalent phase 
functions $ and $' in ()1.4|) will have different values when evaluated for generic g and k. 
However, certain linear combinations of phase functions are independent of gauge. Immedi- 
ately from (jl.4j) . we see that if g leaves k invariant, then $c,(k) is the same in any gauge; 
we call such a quantity a gauge invariant of the first kind. It follows from (jl.2j) that if a 
gauge invariant of the first kind is non-zero (always modulo the integers), then p(k) = 0, 
so there is an extinction in diffraction. For example, let G = 2mm = {e,r,m,rm}, where 
e is the identity, r denotes a 180° rotation, and m is reflection in the x-axis; let L be the 
lattice generated by vectors a and b along the x- and y-axes. Then it is not hard to check 
that $rn(^a + jb) = $rm(^a + jb) = i\, ^e{^) = '^'r(k) = Satisfies the group-compatibility 
condition ()1.3|1 . Thus $m(o) = | is an invariant of the first kind, and the point a is extinct 
in diffraction.'' For further discussion of this example, see §3.21 below. 

This is not the only kind of invariant. Of course any integral linear combination (such as 
$gj(ki) + $g2(k2) where ^fjkj = kj) of gauge invariants of the first kind is still an invariant, but 
for two of the 157 periodic non-symmorphic space groups in three dimensions, the simplest 
gauge-invariant quantity one can construct takes the form 

(2.1) $,(k,)-$,(k,) , 

where g and h commute and where neither term alone is gauge invariant. Konig and Mermin 
|14t IT^ IT^ define the lattice vectors kg and kh in terms of a point q not in the reciprocal 
lattice but with the property that 

(2.2) kg = q — (7q and k/j = q — /iq 

are. The group operations g and h are then elements of the little group of q. We refer the 
reader to their papers or to [2n] for the proof that if the invariant (j2.2j) is non-zero, any 
electronic energy level at wavevector q must be at least two-fold degenerate.^ 

^The point a may also be extinct for the point group Amm on the square lattice (and for analogous star 
lattices 8, 16, . . . ). That there is no systematic extinction for 6mm on the triangular lattice, as a consequence 
of H1.3|l . demonstrates that the calculation of space groups is non-trivial. 

^Equation (|2.1|) is a special case of what we refer to in jH |H| as an invariant of the second kind, which 
always leads to the electronic degeneracy; there also exists an invariant of the third kind, but wc shall not 
need it here. 
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Figure 1. The manifold M, containing two holes, illustrates a topological 
space. AB, C, D, E, F, G, and H are oriented 1-simplices; of these, all but 
AB are 1-cycles. The 1-cycles C and E are 1-boundaries. In particular, C is 
the boundary of sub-manifold S. 

2.2. Topological Invariants^. Figure 1 offers a tour of some of the objects of topology. 
The shaded area represents a two-dimensional manifold M. The region S" is a submanifold; 
the other letters label various oriented curves. An oriented curve is called a 1-simplex. We 
have an intuitive idea, which we shall make formal shortly, of what it means for objects to 
bound one another. As examples, the 1-simplex AB is bounded by the points A and B, but 
the other labeled 1-simplices, which are closed, have no boundary points; such 1-simplices 
with no boundaries will be examples of 1-cycles. Cycles themselves may or may not bound 
submanifolds: the 1-cycle C bounds 5*, but D is not the boundary of any submanifold of M. 
To describe the boundary of M, we must exclude the two holes and so need a combination 
of cycles, which we shall see can be written as the sum H — F — G. Such a formal sum of 
1-simplices is called a 1-chain; any sum, such as C + C + AB = 2G + AB, may be performed, 
although the result might not be very interesting. 

The cycle F encloses the left hole once, while cycles like D and G each enclose the right 
hole once. A cycle might enclose one or both holes an arbitrary number of times; we could 
draw a lemniscate (figure-eight) enclosing the left hole —1 times and the right hole +1 times. 

Boundaries offer a way to identify cycles, such as D and G, that enclose holes the same 
number of times. Since D — G bounds the region between the two cycles, we shall write 
D = G + {D — G) = G, identifying the two cycles because they differ by the boundary 
D — G. Two cycles differing by a boundary will be called homologous; cycles E and G are 
also homologous to each other but not to D and G. A set of homologous cycles forms a 



'The elementary treatment of topology in this section draws on Alexandroff 's slim introduction 
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Figure 2. Figure 2a illustrates the torus T (doughnut) embedded in Eu- 
clidean 3-space. Only the surface forms T. Figure 2b illustrates the same 
torus: the sides with double arrows are to be joined, generating a cylinder, 
and then the sides with single arrows. The triangulation consists of two 2- 
simplices, A and F, three 1-simplices, b, c, and d, and just a single 0-simplex 
(vertex), A. The 1-cycles b and c are also shown in Figure 2a. 

homology class; the enumeration of classes of cycles identified by homology will have a direct 
analogue in the classification of phase functions identified by gauge transformation. 

For this example, each homology class is labeled by two integers: one describing how many 
times each cycle in the class goes around the left hole and another describing how many times 
each cycle in the class goes around the right. Since cycles form an additive group, so do 
homology classes; the homology group of the figure'' is Z x Z. 

The torus, T, of Figure 2a (the surface of a doughnut, not the dough) cannot be deformed 
into the manifold of Figure 1, yet the two have the same homology group. This will follow 
from the fact that a cycle on the torus may enclose the hole any integer number of times 
and that it may enclose the dough (which is also a kind of hole) any integer number of 
times. Since this assertion is perhaps a little less transparent than the corresponding one for 
Figure 1, we start by describing a general method for calculating homology groups. 

First we flatten the torus into the rectangle of Figure 2b. The arrows indicate that opposite 
sides are to be identified, or sealed together: after identifying the two sides with double 
arrows, one obtains a cylinder; identifying the other two sides then gives the torus. Figure 
2b shows the torus with a diagonal added. In this triangulation, the torus T consists of a 
single 0-simplex, or point, labeled A; the three 1-simplices b, c, and d; and the two 2-simplices 
(triangles) A and F. The straight arrows on the 1-simplices indicate their orientations; 2- 
simplices can be oriented as well, as indicated by the curly arrows. The orientations of the 
2-simplices and their constituent 1-simplex sides bear no necessary relation. Neither the 
triangulation of T nor the assignment of orientations is unique; we refer the reader to PP for 
theorems guaranteeing that the homology group is independent of these choices. 

"^The homology group of these curves is called the first homology group; elsewhere ^| we have had 
occasion to employ higher homology groups. We use the common mathematical notations Z for the integers, 
R for the real numbers. 
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Figure 3. The boundary of a 2-simplex (triangle) is the sum of its sides. 
According to ()2.3|1 . we have (xq, xi) — (xq, X2) + (xi, X2) = — (xi, xq) + (xi, X2) + 
(x2,xo), in agreement with the orientations in the figure. 



Labehng a 1-simplex by its endpoints (which may be the same point) and a 2-simplex 
(triangle) by its vertices,^ we define the boundary operator d by 

9(Xo,Xi) = Xi - Xo 

and (9(xo,xi,X2) = (xo,xi) - (xo,X2) + (xi,X2) . 
For example, in Figure 2a, db = dc = dd = 0, and 

dT = -b-c + d 

^^■^^ dA = b + c-d . 

The definitions extend to 1-chains and 2-chains by linearity, and they agree with the 
intuitive notion of boundary. So the 1-simplex (xq, xi) is a 1-cycle — that is, 9(xo, xi) = — if 
and only if xq = Xi. The boundary of a triangle (Figure 3) (xq, xi, X2) is the sum of its 
sides, if we accept the convention that (xq, X2) = — (x2,xo). 

It follows from these definitions that the boundary of the boundary of a 2-chain vanishes. 
Since we may define the boundary of a point (a 0-simplex) as 0, the boundary of the boundary 
of a 1-chain also vanishes. In fact, in the general theory, dd is identically zero. 

The 1-chains whose boundaries are form the group Zi of 1-cycles, while the 1-cycles that 
can themselves be written as boundaries of 2-chains constitute the group Bi of 1-boundaries. 
Homology identifies 1-cycles that differ only by 1-boundaries: that is, we form the quotient 

(2.5) H, = Zi/B^ , 

called the first homology group. 

Applying these ideas to the triangulation of the torus, we first calculate the group Ci of 
all 1-chains. Since there are exactly three 1-simplices, an arbitrary 1-chain takes the form 
ib + mc + nd, with i,m,n G Z; thus Ci is a three-dimensional "space" (properly, module) 



This labeling is not meant to suggest that a simplex is determined by its vertices. The order of the 
vertices does determine the orientation of the simplex. 
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over the integers. In this example, all three 1-simplices have vanishing boundaries, so the 
group of 1-cycles, Zi, is the same as Ci. 

(2.6) Zi = {b,c,d) , 

where the angle brackets indicate the generators (over the integers). From ()2.4|) . we see 
that the only 1-boundaries are integral multiples of b + c — d, so that Bi = (b + c — d), a 
one-dimensional module. As a group, first homology with coefficients in the integers is given 
by 

(2.7) Hi{T, Z) = Z^/Z = ZxZ . 

That is, any 1-cycle is homologous to a linear combination of b (going once around the hole) 
and c (going once around the dough), as we claimed above. 

We have taken the coefficients of simplices to be integers, but we can use any additive 
group. One natural choice is the integers modulo two, Z2 = {0, 1} with the addition rules 
+ = 1 + 1 = 0, + 1 = 1 + = 1. Over Z2, twice any cycle vanishes. This corresponds to 
unoriented simplices: for example, if traversing the 1-simplex C in Figure 1 once is identified 
with traversing the oppositely-oriented simplex once, we have —C = —C + 2C = C, so the 
coefficients are in Z2. If we use unoriented simplices, the homology groups of M and T are 
both two copies of the integers modulo 2: 

(2.8) Hi{M, Z2) = H^{T, Z2) = Z2 X Z2 . 

Before leaving the subject of manifolds, we consider the projective plane, MP^, which can 
be represented, like the torus, as a rectangle, but now opposite sides are identified in opposite 
directions, as in Figure 4.^ The triangulation of the figure comprises the two 0-simplices P 
and Q, three 1-simplices, /, g, and /i, and two 2-simplices, T and ^ . It is easy to see that 
all cycles are generated by h and f + g. We determine the boundaries by 

dT = h-f -g 

^^■^^ d^ = -h-f-g . 

Thus 2h = dT — is a boundary, but the cycle h is not. The other generating cycle, f + g, 
differs from h by the boundary dT; so there exists only a single cycle, which we may take 
to be h, modulo boundaries. The first homology group is 

(2.10) Hi{RP^, Z) = Hi{RP^, Z2) = Z2 . 

The fact that any triangulation can be used to calculate the homology of a manifold is 
related to Euler's formula: 

(2.11) V - E + F = 2, 

where V, E, and F denote the number of vertices, edges, and faces (or 0-, 1-, and 2-simplices) 
in a polyhedron. A polyhedron can be thought of as a polygonal dissection of the sphere; 
if we want, we can triangulate the faces and so obtain a triangulation of the sphere. The 
triangulation (or the polygonal dissection, if one takes the effort to define boundaries of faces 

'^Equivalently, it may be thought of as a disk with opposite points of the circumference identified or, via 
stereographic projection, as a Euchdean plane plus a line at infinity. Since each line contains a point at 
infinity, projection from the focal point P to a line not passing through P is defined on the entire plane 
except for P, hence the term projective plane. 
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Figure 4. The projective plane RP^ (compare Figure 2b). Now, opposite 
sides are to be twisted before being glued together; we despair of showing 
what the result might look like in four dimensions. This triangulation consists 
of two 2-simplices, \E' and T, three 1-simpliccs, /, g, and h, and two 0-simplices, 
P and Q. Twice /i is a boundary, and modulo boundaries, h is the only cycle. 



with more than three sides) can be used to compute the homology groups of the sphere 5, 
and in particular the alternating sum 

(2.12) x{S) = la^nk Ho{S) -rank Hi{S) +riinkH2{S), 

known as the Euler Characteristic of the sphere. It follows from elementary linear algebra 
that the alternating sum is the same whether one takes ranks of homology groups or of chain 
groups: that is, x{S) = V — E -\- F \s the same for any dissection of the sphere. Prom 
the triangulations above, it follows that the Euler characteristics of the torus and projective 
plane are and 1. 



2.3. From cycles to vector fields and de Rham cohomology. Homology of a manifold 
is closely connected to the question of which vector fields are conservative. Let F = Px + Qy 
be a vector field in the plane. If F is the gradient of a potential 0, then the rules of vector 
calculus imply that V x F = and that the line integral of F over any contour vanishes. On 
the other hand, a curlless vector field may or may not be the gradient of a potential. 
The standard counterexample is 

xy - i/x 

It is easy to verify that VxF = but that integrating counterclockwise around the unit 
circle C gives ■ dr — 2n. We could say that F has a 5-function curl at the origin or, 
equivalently, we could cut the origin out of the plane, yielding a "punctured plane," M. In 
the language of homology, the circle C is a cycle in M but not a boundary. If, instead of C, 
we consider the boundary P of a region R that does not contain the origin, then Green's 
Theorem implies that F • dr = 0. 

More generally, consider a vector field that is curlless on a plane punctured any number 
of times. Two cycles that differ by a boundary yield the same line integral, so homology 
naturally describes the classes of line integrals of a curlless vector field. 



(2.13) 
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homology 


cohomology 


l-(co) chain 
l-(co)cycle 
l-(co)boundary 


sum of curves on M 
closed curves (contours for §) 
trivial: bounds region of M 


vector field, F 
V X F = on M 
trivial: §¥ ■ dr = Q 



Table 1. Summary of some of the relations between homology and de Rham cohomology. 



It is equally natural to fix a cycle and to consider different vector fields. Adding the 
gradient of a potential to F does not change any closed line integral, so we identify two 
vector fields that differ by a conservative field. The resulting vector space is the first de Rham 
cohomology group of the punctured plane, M, denoted H^^{M). 

The condition (that F have no curl on M) and the identification (of vector fields differing 
by conservative fields) mirror the cycle and boundary conditions on contours, providing a 
natural duality between de Rham cohomology and Hi{M, M), the homology of M with real 
(not integral) coefficients. On the level of chains and vector fields, the duality is defined by 
the circulation integral, (c, F) = F ■ dr. To see that this is well defined on (co)homology, 
it is enough to notice that ^ F ■ = if c is any cycle and F is conservative (fundamental 
theorem of calculus) or if c is a boundary and F is irrotational (Green's Theorem). 

The language of differential forms |3S1 connects the gradient and the curl. In this language, 
u = F ■ dr = P dx + Q dy is a differential 1-form. The differential of u, giving the curl, 
duj = {dQ / dx — dQ / dy) dxdy, is a 2-form. The exterior derivative d serves as the coboundary 
operator in this context; the vanishing of duo makes uo a 1-cocyde. This resembles the cycle 
condition in homology, except that where d demoted an n-chain to an n — 1-chain, here d 
promotes a 1-form into a 2-form. Similarly, a 0-form is a scalar-valued function, 0, on the 
plane, and d(j) = {dip/ dx)dx+{d(f)/ dy)dy is a 1-coboundary. De Rham cohomology is the space 
of cocycles modulo coboundaries, i.e., vector fields curlless on M modulo conservative fields. ^° 
Table 1 summarizes some of the relations between homology and de Rham cohomology. 

2.4. Crystallographic Invariants as Homology Elements. In the foregoing, we identi- 
fied closed curves (cycles) differing only by boundaries. We then enumerated all the possible 
cycles, up to boundaries. We now describe a similar way to find all crystallographic invari- 
ants, demonstrating the method in the next section. 

First, instead of thinking of a collection of functions ^g, one for each element of the point 
group, think of the phase function $ as acting on a pair consisting of a point-group element 
g and a lattice vector k. We will write this pair as \<-[g]', such an object is an example of 
a 1-chain. The phase function acts on this 1-chain in the obvious way: $(k[5f]) = $g(k). 
Addition is defined on 1-chains as though the various elements of the point group G were 
independent basis vectors in a vector space and with the reciprocal lattice playing the role of 
numeric coefficients.^^ Since the phase function is linear on the lattice, it acts distributively 
over addition. Thus, for example, 

(2.14) ^8L[g]+h[g] + c[h]) = + h)[g] + c[h]) = ^^(a + b) + <I>;.(c) . 



We note that this same example has an equivalent interpretation with contour integrals over functions 
that are analytic in the complex plane except at simple poles. 

^^Since the lattice is not a field, the 1-chains form a module instead of a vector space. 
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Define the boundary operator on a 1-chain by analogy to ()2.3p : 

(2.15) dk[g]=gk-k 
so that for a general 1-chain c = J2i kj 

(2.16) ^c = Y,9i^^-^i ■ 

i 

A 1-chain c is a 1-cycle if dc = 0. Considering the gauge transformation of ()1.4j) . we have 
for the difference between phase functions $ and $' evaluated on the 1-cycle c 

(2.17) <l>'(c)-$(c) = x($^^7k,-k,) = . 

i 

Thus cycles are gauge invariants, and since x is an arbitrary linear function (from the lattice 
to the real numbers modulo unity), any gauge invariant is a cycle. 

Certain gauge invariants are trivial. Applying ()2.15j) . we establish that any 1-chain of the 
form 

(2.18) b = hk[g]~k[gh] + k[h] 

is a 1-cycle. However, any phase function $ evaluated at b must vanish by the group- 
compatibility condition (jl.Hj) . Adding b to any 1-cycle z yields another 1-cycle, z + b, but in 
enumerating gauge invariants, we should not count z and z + b separately. Thus b acts very 
much like the 1-boundary considered in topology. In fact, ()2.18|1 is exactly the definition of 
the boundary of the 2-chain k[5f|/;,]: 

(2.19) dk[g\h] = hk[g]-k[gh]+k[h] . 

Thus, to determine all non-trivial gauge invariants of a group acting on a lattice, we 
calculate the group of 1-cycles and divide by the group of boundaries of 2-chains; the quotient 
group is called Hi{G,L) {G is the point group, L the lattice). As we have noted [22], this 
is less elegant than Mermin's method but, obviating the need for a clever choice of gauge, 
has proven useful in establishing theorems and generalizing results to modulated crystals. 
Moreover, the systematic way in which gauge invariants fall out of this formulation has 
enabled us to find a new type, not previously noted 3j. 

A comparison to the definition of the boundary of the 2-simplex in the second line of ()2.3p 
reveals nearly the same formal structure under the replacement g — >■ {xq,Xi), h (xi,X2), 
gh {xq,X2), the difference being that we considered topological chains over integers, while 
the coefficients in crystallography lie instead in the lattice, on which the group acts non- 
trivially: hence the hk in the first term on the right instead of just k. 

There is a way of visualizing 1-chains in group homology that emphasizes the analogy 
with 1-cycles in topology: picture the 1-chain k[g] as a vector going from k to gk. If g is 
a rotation, it is suggestive to draw the vector along a circular arc, as in Figure 5, in which 
r denotes a 90° rotation. This visualization is consistent with the definition ()2.15|) of the 
boundary of a 1-chain. The formula (j2.19|) for a 1-boundary can be written in the form 
k[gh] = k[h] + (/;,k)[(yf], where = denotes homology as in ^2.21 This is also consistent with 
the vector picture: for example. Figure 5 illustrates the fact that a[r^] = a[r] + (ra)[r^]. 
One shortcoming of this method is that it does not illustrate linearity relations, such as 
(a + b)[(7] = Si[g] + h[g] and {-8i)[g] = -{alg]). 
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Figure 5. It may help to elucidate chains to think of them geometrically; 
a chain k[g] can be pictured as a circular arc pointing from k to gk. The 
boundary 9a[r^] = (9a[r~-^] = a[r] + (ra)[r^], as suggested by the three arcs. 

3. An Example of Calculating Space (Plane) Groups 

To illustrate the homological method, we will calculate the homology group Hi {G, L) 
where the primitive rectangular lattice L is generated by a vector a in the x direction and a 
vector b of different length in the y direction and where the point group G = {e, r, m, rm} is 
generated by a 180° rotation r and the mirror m that leaves a invariant (e is the identity). 
According to the International Tables for Crystallography piJj, the possible plane groups are 
p2mm, p2gg, and p2mg. We shall verify this result in the next subsection; we concentrate 
first on calculating the possible invariants. 

Since our goal is to illustrate a general method as simply as possible, we avoid the shortcuts 
we used in the more complicated example of section 7 of ^H] • In section 14.11 we show how 
the entire calculation can be automated in the symbolic-algebra program, Mathematica jSH], 
for arbitrary examples. 

3.1. Invariants of the Point Group 2mm on the Primitive Rectangular Lattice. 

There are four elements of G, and the lattice is generated by two vectors. To find all the 
cycles, we therefore take boundaries of the eight 1-chains k[(yf] where k is a lattice generator. 
(Any other 1-chain is an integral linear combination of these.) 

db[e\ = 

db[m] = -2b 

2a db[rm] = 

2a db[r] = -2b 



(3.1) 



da[e] 
da [m] 
da[rm] 
dalr] 
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From these eight values, we read off the generators of the additive group of cycles (the 
angle brackets mean "generated by"): 



Zi(G,L) — (zi, Z2, ^3, Z4, ^5, Zq) 



(3.2) 
where 

Zl = b[e] 
(3.3) Z2 = -b[r] + b[m] 

Z3 = b[rm] 

Now, we write down the 32 generators k[(7|/;,] 
2-chains, and calculate their boundaries: 



= a[e] 
= a[m] 
Zq = a[rm] — a[r] 

(k = a, 6 and g, h in the point group) of all 



(3.4) 



bi 


= da[e\e] 


= a[e] 






b2 


= db[e\e] 


= m 








= da[r\m] 


= a[r] 


+ a[m\ 


— a[rm 




= da[m\r] 


= a[r] 


— a[m\ 


— a[rm 


h 


= db[r\m] 


= —b[r] 


+ b[m] 


— 6[rm] 




= db[m\r] 


= b[r] 


— 6[m] 


— b[rm] 



All of the remaining boundaries are linear combinations of these six boundaries. Furthermore, 
it is evident that the six boundaries in (j3.4j) are integrally linearly independent, so^^ 

(3.5) 5i(G,L) = (61, 62,63, fe4,fc5,&6) ■ 

We expect and verify that every boundary in ()3.5p can be written as an integral linear 
combination of the cycles in ()3.3|) (since all boundaries are cycles). However, the cycles zi 
and Z4 are actually boundaries, so we throw them out. Furthermore, — zq = b^, which is 
a boundary, so we write z^ = zq (equality up to boundaries); similarly, Z2 = Z3. This leaves 
only two cycles, z^ and z^, which are obviously linearly independent and not boundaries. 
Finally, we note that 22:3 = — (65 + b^) and 2^5 = 63 — 64 are boundaries. Thus there are only 
four invariants: 

(3.6) Hi{G,L) = {0, Z3, z,, ZS + Z5} . 

3.2. Connection between invariants and space groups (cohomology). Now we're 
ready to classify plane groups. According to the Rokhsar-Wright-Mermin specification, we 
need to find all functions $, linear modulo unity on the lattice, satisfying (|1.3|) . Since 
22:3 and 2z5 are boundaries, we must have 2^{z3) = 2^(25) = (all modulo unity), so that 
there are four possible gauge-invariant values for the phase function, given in Table 2. 

We now make two assertions, the proof of which is the main content of reference j^E]- 
First, two phase functions not related by a gauge differ in their values on at least one gauge 
invariant in Hi{G,L). Thus, there are no more than four space groups in the example. 
Second, for each possibility, a phase function exists. 



'in M.2\ we prove that Zi and Bi always have the same rank. 
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possibility 


$(2:3) = 


$(^5) = $^(a) 


1 








2 


1/2 





3 





1/2 


4 


1/2 


1/2 



Table 2. The four sets of gauge-invariant values a phase function might take 
for the point group 2mm on the primitive-rectangular lattice. These cor- 
respond to the plane groups p2mm, p2mg, p2gm (which is just a different 
setting of p2mg), and p2gg. 





homology 


cohomology 


1- ( CO ) chain 
1- ( CO ) cycle 
l-(co)boundary 


gauge invariant 
trivial by group compatibility 


$ 

satisfies group compatibility 
trivial by gauge transformation 



Table 3. Summary of the application of group homology and cohomology to 
Fourier-space crystallography. (For additional applications and higher-order 
(co)homology, see [26, 9J.) Symbols: kj are reciprocal-lattice vectors and gi 
point-group elements. $ is a phase function. Compare Table 1. 



The proof in [211 i"ests on the duality of homology and cohomology, in which $ plays the 
same role with respect to a chain ^ ]<.[g] as the vector field F played with respect to a contour 
in subsection 12.31 Think of $ as acting on g to give $p; $g acts linearly on k G L to give 
the number $g(k), modulo unity as always. The 1-cochains form a group under addition. A 
1-cochain is a 1-cocycle if its coboundary, defined by 

(3.7) {d^){g,h) = <!>goh-^,h + ^h , 

vanishes: this is simply the group-compatibility condition on $, so we are clearly interested 
in 1-cocycles. Note that where the boundary operator d decreased the number of group 
arguments, the coboundary operator d increases it. A linear function x on the lattice is 
called a 0-cochain; its coboundary is given by 

(3.8) {dx){9) = X° 9 -X ■ 

Two different 1-cochains $ and $' that differ by a coboundary of the form of ()3.8|) are related 
by a gauge function, according to ()1.4j) . so we identify them. The group of 1-cocycles modulo 
1-coboundaries, or first cohomology, is labeled H^{G, L), where L denotes the group of linear 
maps from L to M/Z; if x is an element of L, then 27rzx(k) is a choice of phase for each k 
among a generating set for the lattice L. Table 3 summarizes homology and cohomology; 
note the similarity between the requirement for homology and triviality for cohomology, and 
vice-versa. 

The task of classifying space groups, then, comes down to calculating the cohomology 
group of cocycles modulo coboundaries, and this group is dual to the homology group of 
cycles modulo boundaries. Thus, the four possibilities in Table 2 above are exactly the four 
possible space groups for this example. 
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However, there are only three space groups for 2mm on the primitive rectangular lattice: 
p2mm, p2mg, and p2gg. Possibilities 2 and 3 in the table simply exchange the x and y 
directions. 

3.3. Comparison to torus and lemniscate. Table 2 shows that the homology (or co- 
homology) group of the point group 2mm acting on the primitive-rectangular lattice is 
isomorphic to Z2 x Z2, the same homology group we considered in two examples in section 
12. 2[ ^^ The topological examples contained two holes in the sense of closed curves that did 
not bound: Figure 1 without orientation has two literal holes, while the doughnut has the 
hole in the center and the fried dough (which is not part of the surface and so is no different 
from a hole). In each case, exchanging the two holes reduces the number of combinations 
from four to three. The close analogies tell us that we can consider crystallography, as well, 
in terms of topological spaces. A point group and lattice admitting no non-symmorphic 
space groups is a trivial space, in which every closed curve can be collapsed to a point, while 
non-symmorphic space groups are non-zero cohomology classes in spaces with one or more 
holes. We shall show in section 14.11 that in every case, the invariants correspond to some 
number of holes, each of which can be lassoed by a cycle a finite number of times before it 
vanishes. (More prosaically, the (co)homology group is isomorphic to the direct product of 
some number of cyclic groups of varying order.) 



4. Generalizing the calculation 

4.1. Automation of calculation in Mathematica. The algorithm of section 13. H re- 
quiring no clever choices of gauge or anything specific to the point group or lattice, is 
easily automated. Our Mathematica implementation, available at the preprint archive, 
"cond-mat/0301601 , takes only 13 lines of substantive code. A quick tour of the algorithm 
may assist those trying to do the calculation by hand as well as users of computer-algebra 
packages. 

(1) The point-group generators and their actions (to the right) on the lattice are specified 
by square matrices of dimension equal to the rank of the lattice. For example, the 
two-fold rotation r and mirror m in section 13.11 are the matrices 

"=(~J -?)' "^=(J -1) 

acting on the column vectors a = (q) and b = (°). 

(2) The group Ci of all 1-chains is generated by the combinations \<-[g], where k is a 
lattice generator and g E G. The computer stores all these combinations in the 
list c. Think of this list, say c = (ci, . . . , c,,), as a column vector. 

(3) Take the boundary (j2.16j) of each Cj. The NullSpace [ ] function in Mathematica, 
applied to the list of boundaries, gives a matrix z that expresses generators Zi of the 



^''Once cohomology has been calculated, it is still necessary to consider rotations and (in quasiperiodic 
cases) scale invariance of the lattice. See |H1 118| for a full accounting of Bravais class, arithmetic crystal 
class, and space-group type. 

^^To the manifold with two holes and the torus, each with coefficients in Z2, we may add RP^ x RP^ 
with coefficients in Z. 



(4.1) 
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cycle group in terms of c: 



(4.2) 



(4) The group of all 2-chains is generated by the combinations k[(7|/;,], where k is a 
lattice generator and g, h E G. The boundaries hi of these 2-chains ()2.19j] generate 
the group Bi of 1-boundaries. Express these generators in terms of c: 



(4.3) 











1 = bz3 j 











(5) Find a rational transformation z2 that transforms a cycle from the basis of Ci into 
the basis of Zi. Our package uses Pseudolnverse [ ] . 

(6) From steps El and El the product bz3.z2 expresses the generators hi of Bi in the basis 
of Zl. We call LatticeReduce [ ] to find the smallest set of rows that will generate 
Bi over the integers, yielding the matrix b. Thus 



(4.4) 



gives a minimal set of generators for Bi. We explain in 14.21 below that Bi and Zi 
have the same rank, so that b is a square matrix. We can calculate the first homology 
group Hi = Zi/Bi as follows. The Smith normal form Thm. 25.26], |12j of the 
matrix b is a diagonal matrix d such that there are invertible integral transformations 
p and q with 

(4.5) d = pbq . 

Thus q and p describe new bases for Zi and Bi, say {z^} and {&•}. In these new 
bases, h'^ = di z[\ in other words, rfj times cycle z[ is a boundary. Thus the homology 
group B\ can be described as all linear combinations of z^, . . ., z^, where the coefficient 
of z[ lies between and di — 1, and there are no relations other than diz\ = 0. 
(7) Our function SpacegroupHlFull [ ] returns the matrices c, z, b, d, p, and q above. 

Returning again to the example of sectioning we find 

(4.6) d = diag(l,l,l,l,2,2) . 

The module Zi has dimension six, as in ()3.3|) . The four unit entries indicate that the 
corresponding cycles are boundaries, while the two entries of 2 tell us that two, z'^ and z'q, 
are not, although 2zi, and 2zq are. We invert the Smith transformation to solve for z'^ and 

4-- 

(4.7) p-Miag(0,0,0,0,l,l) q"i . 

This yields a matrix with two non-zero rows, (0,-1,0,0,0,0) and (0,0,0,0,0,1), giving 
zi, = —Z2 and Zq = zq. The z matrix converts the Zi basis into the basis of Ci, which we 
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(* For g a list of matrices, return <g>, the group generated. *) 
groupgen[g_] := FixedPoint [Union [Flatten [Outer [Dot ,#,g, 1] , 1] ,#] &, g] 

(* Project a term of form (vl y[gl]) onto elements of cl. *) 
clproj[x_,cl_] := ( (x/ . {y [g_] ->If [g==# [ [2] ] , 1 ,0] }) .# [[1] ] )& /@ cl 
latnum = Dimensions [x [ [1] ][ [1] ]][ [1] ] ; (* x = list of generators *) 
latgen = IdentityMatrix [latnum] ; 
group = groupgen[x] ; 
12] c = Distribute [{latgen, group}, List]; 

Gl z = NullSpace [Transpose [(# [[2] ] .#[[1]] - #[[!]])& /@ c]]; 
131 c2 = Distribute [{latgen, group, group}, List]; 

bz2 = Union[((#[[3]] .#[[l]])y[#[[2]]] - # [ [1] ] y [# [ [2] ] .# [ [3] ] ] + 

#[[l]]y[#[[3]]])& /@ c2] ; 
bz3 = Union [ (clpro j [#,c] )& /@ bz2] ; 
121 z2 = Chop [Rationalize [Pseudolnverse [N [z] ]]] ; 
1^ b = LatticeReduce [Rationalize [Chop [bz3 . z2] ]] ; 
{d,{p,q}}=ExtendedSmithForm[b] ; 
Table 4. The substantive lines of Mathematica code described in section I^^Tl 
The circled numbers refer to the steps listed in the text. The full version of 
the code is available at the preprint archive, |cond-mat/030160l] 



decode with the help of the list c to compare with the results of section ITTl 

(4.8) 1 : : 

We admit this algorithm to be inelegant: like direct-space formulations, it requires 
the construction of matrices of dimension equal to the rank of the lattice. Unlike the 
Rokhsar-Wright-Mermin method, it cannot treat related point groups, e.g., 2^ mm = 
{4mm, 8mm, (16)mm . . . }, all at once but rather requires a separate calculation for 
each. However, we find it useful in those cases (such as 46-fold symmetry and modulated 
lattices) where discovering a good choice of gauge, requisite to Rokhsar-Wright-Mermin, 
seems too difficult. 

4.2. Finite number of Space Groups. When applying the Rokhsar-Wright-Mermin ap- 
proach to quasicrystals, there are a few surprises. One is that, even in two dimensions, a 
finitely-generated lattice that is not discrete may have an infinite point group 35 . Assume, 
as we have been doing implicitly so far, that the point group is finite. Happily, this is enough 
to guarantee that there are only finitely many space groups associated to the point group G 
and the lattice L. As explained in ^3.21 this amounts to saying that the homology group 
Hi{G, L) is finite. In any example, this can be checked by verifying that the matrices b and d 
of ^4.11 have the same rank as the cycle group Zi. 

In order to prove generally that Hi{G, L) is finite, it suffices to show that Nz is a boundary 
for any 1-cycle z, where is the order of the point group G. Grant this for the moment. 
There are then finitely many 1-cycles, say zi, Z2, . ■ ■ , z^, (where M is the rank of the cycle 
group Zi) that generate all the others: any 1-cycle z can be written as a linear combination 
of the generators Zi. In order to represent all homology classes, it suffices to take coefficients 
between 1 and A^, since Nzi is a boundary. Thus there are only finitely many homology 
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classes. We can also explain this in terms of matrices: if Nzi is a boundary for each gen- 
erator Zj, then b'b = NIm for some integer matrix b', where Is denotes the s x s identity 
matrix. This implies that b has rank s. 

It remains to explain why times a cycle is a boundary. Let z = Ylh^h[h] be a cycle, 
where h runs through the point group. According to ()2.16p . the cycle condition dz = can 
be written in the form 

(4.9) J2hkn = Y,^h. 

h h 

Construct a 2-chain c such that dc = Nz using the fundamental technique of averaging over 
the group: let 

(4.10) c = J2M9H 
According to ^TH^ . 

(4.11) dc = J^ih^hM -J2^h[9h] + J2^H[h]. 

g,h g,h g.h 

The first two terms in (j4.1H) cancel: to see this, fix g and apply the cycle condition (j4.9j) in 
the first term; and fix h and replace g with gh~^ in the second term. The summand in the 
third term in (j4.1H) is independent of 5^, so the third term is simply Nz. 

4.3. From Real Space to Fourier Space. So far, we have described the Fourier-space ap- 
proach to crystallography in the language of (co)homology, explained how this language lends 
itself to explicit calculations, and pointed out the analogy with topological (co)homology. 
The real advantage of adopting this framework comes from applying the theorems developed 
over the years for homological algebra in general and group cohomology in particular. The 
duality between phase functions and crystallographic invariants f ^2.4|) is a special case of 
one such theorem, the duality between homology and cohomology. The calculation of invari- 
ants {^^ 13. II and 14. 1|) is one example of a standard technique. The finiteness result f ^4.2|) is 
another. 

In this subsection, we summarize the connection between the real-space and Fourier-space 
approaches to crystallography, using the common language of group cohomology. Each step 
in this comparison involves a standard result, explained in textbooks such as [S], but we 
will not try to reproduce these explanations. At the end, we illustrate these generalities by 
returning to the example of the point group 2mm. 

In real-space crystallography, one considers a crystal and its space group Q. For a periodic, 
(i-dimensional crystal, the space group consists of all isometrics of that preserve the 
crystal. If we are interested in a quasiperiodic crystal, we construct a periodic crystal in a 
higher-dimensional space, say M^, that projects to the quasiperiodic one; for this discussion, 
Q will be the space group of the higher-dimensional, periodic crystal. (In the periodic case, 
take D = d.) Let T denote the lattice of pure real-space translations in Q. 

In group-theoretic terms, the point group is the quotient G = Q /T . It can be thought of 
as the group of "macroscopic symmetries" : think of all translations as "microscopic," and so 
identify any two elements of the space group if they differ by a translation. Both G and T 



Different choices of the periodic ZJ-dimcnsional crystal may lead to different ZJ-diniensional space groups 

m 
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are fairly easy to describe: G is a finite subgroup of the orthogonal group and T consists 
of all integral linear combinations of D generating vectors. These two groups are also fairly 
easy to determine experimentally. It remains to describe the space group Q in terms of 
these two simpler groups. One possibility is that the point group is contained in and the 
space group is generated by G and T. In this case, the space group is called symmorphic in 
crystallographic terminology, or a semidirect product in the language of group theory. 

Not all space groups are symmorphic. A standard result of group theory states that the 
space groups corresponding to G and T are classified by the cohomology group^^ H'^{G,T). 
Using the Long Exact Sequence of cohomology (c/. the last Remark in §5 of [H]), this is 
isomorphic to H\G,R^/T). 

In the periodic case, the Fourier lattice L is the dual of the translation group T. That 
is, the Fourier lattice consists of those vectors k such that k • t is integral for all t in 
the direct lattice T. (Those who worry about where the factor of 27r went should refer 
to and ()1.2p .) Drager and Mermin |Hj realized that this can be turned around in 

the quasiperiodic case: the translation group T and the super-space can be described, 
without explicit coordinates, as duals of the Fourier lattice. Algebraically, this leads to the 
isomorphism MP /T = L, where L denotes the linear maps from L to M/Z as in ^3.21 Thus 
space groups are classified by the cohomology group H^{G,L), which is exactly the set of 
phase functions (identifying those that differ by a gauge transformation) described in S ^l.ll 
and EH 

We now return to the point group 2mm as in ^3.11 Let T denote the lattice in dual 
to the Fourier lattice L generated by a and b. We will use the terminology of cocycles and 
coboundaries introduced in ^3.2( with $ replaced by t. We will also use the notation {g\t} 
for elements of the space group, where g is in the point group and t is a translation (not 
necessarily in T): this acts on the vector x by the formula {^fltlx = (yfx + 1. It follows that 

{fi'ilti} ■ {fi'2|t2} = {9l92\tl + fi'i(t2)}. 

For each element g oi G = {e, r, m, rm}, choose a translation tg such that {(/Itg} is in the 
space group. By choosing the origin appropriately, we can arrange it so that hes along the 
x-axis and trm lies along the y-axis. For each of the glide reflections {m|tm} and {rm|trm}, 
there are two possibilities: the translational part can lie in T, or it can be half of a translation 
in T. Since m and rm generate the point group, the choice of and t^m determines the 
space group. We shall see that the different possible choices lead to the four possibilities in 
Table 2. We shall postpone the veriflcation that each of the four choices actually leads to a 
space group with translation subgroup T and point group G. 

The definition of multiplication in the space group shows that the combination (c/. fj3.7|) ) 

(4.12) dt{g,h)=tg-tgj, + g{tn) 

is the translational part of {g\tg}{h\th\{gh\tgii}'^ and so lies in T . In fact, dt is a 2-cocycle 
with values in T }^ The choice of origin, mentioned above, does not change the value of dt. 
Choosing a different translation tg for each g would change dt by a 1-coboundary, so it 



"The real-space approach cannot be used to describe quasicrystals with infinite point groups, but it is 
not clear that these are physically interesting. 

^^This result applies because T is a normal, Abelian subgroup of Q. 

^^li we allow cochains with values in K^, then tg is a 1-cochain and dt is its boundary. It is therefore 
automatically a cycle. However, dt is not necessarily the boundary of a 1-cochain with values in T. 
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makes sense to identify cocycles that differ by a coboundary. Thus dt determines a class in 
H'^{G,T), and this cohomology class corresponds to the given space group. 

Another point of view is that one should not choose a particular but should consider 
all possible choices at once. In other words, instead of the vector tg, consider its coset 
tg = tg + T, as an element of the quotient group R^/T. From this point of view, t is a 
1-cochain with values in R^/T. The fact that dt takes values in T means that t is a cocycle. 
A different choice of origin changes this cocycle by a coboundary. Thus the space group is 
described by the class of t in R^/T). The relation between t and dt illustrates the 

general isomorphism between H^{G,M.^ /T) and H^{G,T). 

The phase function that describes the space group is simply $g(k) = k • t^, where the dot 
product defines the duality pairing between Fourier and real space. Note that changing 
by a translation in T changes the phase by an integer, which we always ignore. This illus- 
trates the general isomorphism between H^{G,MP /T) and H^{G, L). As promised, the four 
possible choices of and t^^n indeed correspond to the four possibilities in Table 2. From 
either the Fourier-space point of view or the real-space one, two of the four space groups are 
equivalent under exchanging and rescaling the axes. 

We already know, from ^3.21 that H^{G, L) contains exactly four elements. Therefore, the 
isomorphism between H^{G, T) (which classifies space groups from the direct-space point of 
view) and H^{G,L) (the group of phase functions up to gauge equivalence) gives a round- 
about verification of the point we omitted above: all four ways of choosing the glide reflections 
tm and trm actually lead to space groups with translation subgroup T and point group G. 

4.4. Is there a topological space for every crystal class? We have offered an anal- 
ogy between the classification of space groups^^ by group cohomology and some topological 
problems on manifolds. One can ask how deep the connection really goes. We might wish 
to know whether, for a point group G with its action on a lattice L, we could always find a 
topological manifold M and coefficient group C such that 

(4.13) H,,{G,L)^H^{M,C) 

for all integers n > 0. We do not know, in general, whether such a manifold exists. However, 
there does exist a general construction for a topological space so that (j4.13p works for n = 1. 
Since Hi{G,L) is isomorphic to a finite direct product of cyclic groups — -ffx(Z;^,Z), 
where denotes the lens space jTl §24.4], the corresponding direct product of lens spaces 
has the desired first homology group Prop. 0.8]. 

The topological view provides a picture of what one is doing in the Rokhsar-Wright- 
Mermin formulation. There, the task of finding all gauge invariants is complicated first by 
the possibility that these gauge invariants may be linear combinations of phase functions, as 
in ()2.H) . and second by the need to find a clever choice of gauge to make all values that are 
not invariant, or derivable from them, vanish. Viewed cohomo logically, the second point is 
no concern, since the gauge-invariant values completely determine the phase function 
The first complication is also no concern in the homological approach; the fact that (j2.ip is 
a linear combination of terms simply reflects an inopportune choice of basis vectors for the 
cycle group Zi. A more convenient basis would include an invariant like \<.h[g] — kg[/i], and 
each basis element can be thought of as corresponding to a hole in a topological space. 



That is, just that part of the classification that concerns famihes of gauge-equivalent phase functions 
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Figure Captions 



Figure 1. The manifold M, containing two lioles, illustrates a topological space. AB, C, 
D, E, F, G, and H are oriented 1-simplices; of these, all but AB are 1-cycles. The 1-cycles 
C and E are 1-boundaries. In particular, C is the boundary of sub-manifold S. 

Figure 2. Figure 2a illustrates the torus T (doughnut) embedded in Euclidean 3-space. 
Only the surface forms T. Figure 2b illustrates the same torus: the sides with double 
arrows are to be joined, generating a cylinder, and then the sides with single arrows. The 
triangulation consists of two 2-simphces, A and F, three 1-simplices, b, c, and d, and just a 
single 0-simplex (vertex), A. The 1-cycles b and c are also shown in Figure 2a. 

Figure 3. The boundary of a 2-simplex (triangle) is the sum of its sides. According to 
()2.3p . we have (xq, Xi) — (xq, X2) + (xi, X2) = — (xi, Xq) + (xi, X2) + (x2, Xq), in agreement with 
the orientations in the figure. 

Figure 4. The projective plane MP^ (compare Figure 2b). Now, opposite sides are to be 
twisted before being glued together; we despair of showing what the result might look like in 
four dimensions. This triangulation consists of two 2-simplices, and T, three 1-simplices, 
/, g, and h, and two 0-simplices, P and Q. Twice h is a boundary, and modulo boundaries, 
h is the only cycle. 

Figure 5. It may help to elucidate chains to think of them geometrically; a chain k[g] 
can be pictured as a circular arc pointing from k to gk. The boundary (9a[r^] = c}a[r^^] = 
a[r] + (ra)[r^], as suggested by the three arcs. 




Figure 2 a. 




Figure 2b. 




Figure 3. 
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Figure 5. 
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